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ABSTRACT 


We study the distortions induced by peculiar velocities on the redshift-space 
correlation function of galaxies of different morphological types in the Pisces-Perseus 
redshift survey. Redshift-space distortions affect early- and late-type galaxies in 
different ways. In particular, at small separations, the dominant effect comes from 
virialized cluster cores, where ellipticals are the dominant population. The net result is 
that a meaningful comparison of the clustering strength of different morphological types 
can be performed only in real space, i.e., after projecting out the redshift distortions 
on the two-point correlation function £(r p , it). A power-law fit to the projected 
function w p (r p ) on scales smaller than 10 h~ l 2 3 * 5 Mpc gives ro = 8.35lg;7g h _1 Mpc, 

7 = 2.05lo'()g f° r the early-type population, and ro = 5.55 +q 45 h -1 Mpc, 7 = 1.73+o;gg 
for spirals and irregulars. These values are derived for a sample luminosity limited to 
Mzw < —19.5. We detect a 25% increase of ro with luminosity for all types combined, 
from Mzw = —19 to —20. In the framework of a simple stable-clustering model for 
the mean streaming of pairs, we estimate <xi2(l), the one-dimensional pairwise velocity 
dispersion between 0 and l/i -1 Mpc, to be 8651^5 km s _1 for early-type galaxies and 
345 tH kms 1 for late types. This latter value should be a fair estimate of the pairwise 
dispersion for “field” galaxies; it is stable with respect to the presence or absence of 
clusters in the sample, and is consistent with the values found for non-cluster galaxies 
and IRAS galaxies at similar separations. 
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1. INTRODUCTION 


Peculiar velocities distort the maps of the galaxy distribution when redshifts are used as a 
measure of distance through the Hubble relation. The observed distortions contain important 
information on the statistical properties of the large-scale motions of galaxies, presumably due to 
the gravitational influence of the true underlying mass distribution. In particular, the two-point 
correlation function in redshift space £(s) differs from that in real space £(r) in two respects: on 
small scales correlations are suppressed due to the virialized motions in rich clusters, which in 
redshift space elongate structures along the line of sight; on large scales coherent motions produced 
by infall into over dense regions or by outflow out of under dense regions enhance correlations. 

Galaxies of different morphological types inhabit different environments, following a 
well-established morphology-density relation (e.g., Dressier 1980; Postman &; Geller 1984). As 
a consequence, they display significantly different clustering properties (Davis & Geller 1976; 
Giovanelli, Haynes & Chincarini 1986; Iovino et al. 1993; Loveday et al. 1995; Hermit et al. 
1996): ellipticals and SO’s dominate dense cluster cores and are therefore more clustered than 
spirals and irregulars. However, their association with the deep potential wells of clusters implies 
that they have higher peculiar velocities on average, so that the small-scale £(s) for early-type 
galaxies is more strongly suppressed with respect to £(r) than for late types. The consequence 
is that a correct comparison of the clustering properties of different morphological types requires 
understanding in detail their respective redshift space distortions. One way to avoid this problem 
is to measure the angular correlation function w(6) on two-dimensional catalogues (Giovanelli et 
al. 1986; Loveday et al. 19950). The effects of redshift space distortions also need to be quantified 
in comparing the angular clustering of distant objects with the clustering in redshift space of 
galaxies at low redshift (e.g., Iovino et al. 1996). 

The standard way to quantify redshift distortions is to split the separation vector of a pair 
of objects into components lying on the plane of the sky, r p , and along the line of sight, n r, 
and to compute the correlation function ^(r p ,7r) as a function of these two components. The 
iso-correlation contours of £(r p ,7r) will be stretched along the n direction at small separations, 
due to the effect of large velocity dispersions, and compressed at large scales, as a consequence 
of large-scale coherent motions. Projecting £(r p , 7 r) onto the r p axis gives the projected function 
w p ( r p), which is independent of redshift distortions and can be directly expressed as an integral 
over the real-space correlation function £(r). The equation relating w p {r p ) and £(r) thus allows 
one to recover the latter via direct inversion or modelling. 


On the other hand, modelling the distortions of £(r p , n) allows one to characterize the pairwise 
velocity distribution function. We are particularly concerned with the second moment of this 
distribution function, ui 2 {r ), in this paper. Davis & Peebles (1983, hereafter DP83|) used the 


'These latter authors perform also a cross-correlation between morphological subsamples of the APM-Stromlo 
redshift survey and the APM angular catalogue, see § 4.2 







- 3 - 


CfAl survey ( Huchra et al. 1993 ) to measure < 712 ( 1 ), the value of ^(r) at r = 1 /i -1 Mpc, to be 
340 ± 40kms^ 1 . While the analysis of the IRAS 1.2 Jy redshift survey by Fisher et al. (1994b, 
hereafter F94b ) produced a similar result, < 712 ( 1 ) = SWl^gkms^ 1 , re-analyses of the CfAl survey 
( Mo, Jing fe Borner 1993] , |Zurek et al. 1994 , [Somerville, Davis fe Primack 1996 ), and of larger 
optical redshift surveys have shown a large range of values, going as high as lOOOkms -1 (Marzke 


et al. 199q, |Guzzo et al. 199q, [Lin 19951). In particular, < 712 ( 1 ) is found to be quite sensitive to the 


presence or absence of one or two rich clusters, even in volumes as large as the CfA2. It seems 
plausible that while the CfAl value was strongly affected by the smallness of the volume surveyed, 
that measured from the IRAS 1.2 Jy survey reflects the specific nature of IRAS galaxies, which 
are mostly star-forming, late-type galaxies which are under-represented in rich clusters relative to 
optically selected galaxies ( Strauss et al. 199^ ). 


The case of IRAS galaxies explicitly illustrates the dangers of using a specific class of objects 
to draw conclusions on statistics of the velocity field, in particular at small separations: the answer 
depends sensitively on the morphological type of the tracer we are using. This paper addresses 
this issue in detail: 1) What is the difference in the clustering strength of early- and late-type 
galaxies measured in real space? 2) What is the small-scale velocity dispersion for the two classes 
of objects? 

The outline of the paper is as follows: in § |2[ we present the data we will use for our analyses. 
We discuss the measurement of £(r p ,ir) in § and present our results in § [|. Our conclusions are 
summarized in § 


2. THE DATA: DEFINITION OF THE SAMPLES 


We use the Perseus-Pisces redshift catalogue (cf., piovanelli fe Haynes 1991 ), which includes 
redshifts for all Zwicky galaxies QZwicky et al. 1961—1968|) in the positive-declination part of the 


South Galactic Cap (i.e., about 21 ft < a < 5 , 0° < 5 < 50°). As Giovanelli, Haynes fe Chincarini 


1986| make clear, the Perseus-Pisces redshift survey is affected by Galactic extinction around the 
edges. For statistical studies, therefore, it has to be properly restricted. The Zwicky catalogue 
is nominally complete to mz w = 15.7: we thus impose an extinction-corrected magnitude cut of 
15.5, trim the survey to 22^ < a < 4 h , 0° < S < 42°, and apply the additional cut indicated by 
the heavy line of Figure [I]. This excludes nearly all regions with absorption Ab > 0.2, as given by 
the extinction maps of |Burstein fe Heiles (1978) , while it leaves in the core of the Perseus cluster 
(a ~ 3.2^, 5 ~ +41°) to allow us to study the robustness of our results to the presence of the 
richest cluster in the region. The magnitude-limited sample selected in this way contains 4111 
galaxies. One potential problem of our selection criteria is that they push the Zwicky catalogue 
to its completeness limit. In particular, with the extinction correction we include galaxies with 
observed Zwicky magnitudes mz w ~ 15.5, where magnitude errors are large (e.g., Bothun & 
Cornell 1990). We shall show in § 4.2 that our principal results are indeed quite robust to these 
uncertainties in the parent photometric catalogue. 
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All the analyses of this paper are done with volume-limited subsamples of the data. That 
is, we select a lower-limiting luminosity (or equivalently, an upper limit in absolute magnitude), 
and a corresponding maximum distance implied by our apparent magnitude limit, giving us 
uniform sampling throughout the volume. This has the effect of de-emphasizing the Pisces-Perseus 
chain relative to the magnitude-limited case, because in the latter, the selection function peaks 
near the redshift of the supercluster. Also, this choice is crucial for discussing luminosity 
effects, and eliminates uncertainties related to the weighting schemes necessary when analysing 
magnitude-limited samples. 

Table |I] summarizes the parameters of the volume-limited samples we have used. The range 
of absolute magnitudes covered by the subsamples reflects the need to maximize our volume, 
while keeping a sufficient number of objects within it. Absolute magnitudes were calculated 
assuming Hq = 100 kms _1 Mpc _1 . The E-19.5 sample (“ellipticals”) contains galaxies with 
early morphological types (E, SO, and SOa), while the S- samples (“spirals”) contain all galaxies 
classified as spirals or irregulars. The morphological information available in the catalogue is in 
reality quite finer, subdivided into 14 classes. The morphological coding is is from the UGC for 
those galaxies in that catalog, and has been estimated from sky survey plates for the remainder 
of the galaxies (Giovanelli, Haynes & Chincarini 1986). To maximize the statistics within the 
volume-limited samples, we restrict our analysis to the two broad groups of early- and late-tape 
galaxies. The spiral class is however large enough to define several samples to different absolute 
magnitudes (S—19, S—19.5, and S-20); we also define a spiral sample trimmed to exclude the 
Perseus cluster (S-19.5-NOP). Finally, we define an equivalent sample including all morphologies 
for comparison, PP-19.5-NOP. 

The cone diagrams of Figure ^ show the galaxy distribution in PP-19.5 and PP-20, in which 
all morphological types are included, while Figure ^ shows the corresponding distribution of the 
E-19.5 and S-19.5 samples, respectively. Note that the Perseus-Pisces chain, the overdensity at 
cz ~ 6000kms _1 , is more prominent in the ellipticals than the spirals, consistent with the observed 
distribution on the sky ( Piovanelli, Haynes &: Chincarini 1986| ). This is made quantitative in 
Figure |i|, which shows the redshift histogram of each of these subsamples; the dashed line in each 
case, proportional to r 2 , gives the expected distribution in the absence of structure. 


3. ESTIMATING THE TWO-POINT CORRELATION FUNCTION IN 

REDSHIFT SPACE 

3.1. Definitions of r p and n 

The effect of redshift-space distortions can be understood through the correlation function 
£(r p , 7 r), where the radial separation of pairs is split into two components: n, parallel to the line of 
sight, and r p , perpendicular to it. There are two definitions of these quantities in the literature. 
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Given two galaxies at redshifts v\ and v 2 , separated by angle 6, DP83 define: 


7 r 


I 


1 / v 0 

\vi-V 2 \ r p = — {v i + n 2 )tan- 
n 0 2 


(1) 


Note that the quadrature sum of r p and 7r is not equal to the redshift space distance s between the 
galaxies. Recognizing this, Fisher et al. (1994a, hereafter F94a ) use a slightly different definition. 
They define the line of sight vector 1 = (vi + vo)/2 and the redshift difference vector s = vi — V 2 , 
leading to the definitions 

' 2 ( 2 ) 


s • 1 


7T = 


H 0 \l 


r 2 = 
p ~ 


s • s 


H 2 0 


— 7T 


If we recast the F94a formulation in terms of 6, we find: 


7 T = - v 2 \ + 0(d 2 ) 

-Ho 


- ^ ( y i + v 2 ) tan 


Av\V 2 


(3) 


(ui + v 2 ) 2 ’ 

which shows explicitly that the two definitions are not strictly equivalent, even in the small-angle 
approximation. For our analysis we use Eq. (|2j), but checks using Eq. ([!]) show only negligible 
differences in our £(r p ,7r) maps. We conclude therefore that one can make direct comparison of 
results obtained with the two different definitions. 


3.2. Measuring £(r ; 


p; 


vr 


We estimate the quantity ^(r p , 7r) using the method of |DP83| . A catalog of ur = 100,000 
uniformly distributed points with the same boundaries as the real sample is prepared. We count 
the number of pairs in 1 h~ l Mpc bins of r p and 7r among the no galaxies [Ngg^p-, 7 t)], and 
between the galaxies and the random sample [NGR(r p , 7r)]. Our estimate of the correlation function 
is then 


i{r p ,TT) 


N G G{r pi 7r) 2n R 
NGR{r p ,TT) n G 


(4) 


Because our samples are volume-limited, each galaxy gets equal weight, and thus we do not 
apply the statistical weights needed when analyzing magnitude-limited catalogues. For the most 
part, we confine ourselves to scales less than 10 h~ 1 Mpc, and thus there is little benefit to using 
the alternative estimator of Hamilton (1993). This is less sensitive than is Eq. (Q) to uncertainties 
in the mean density, and therefore is important when measuring the correlation function on very 
large scales. 


3.3. Error Estimation and Maximum—Likelihood Fits 


Following Ling et al. (1986) , we use bootstrap resampling to compute statistical errors on 
our estimates of quantities we derive from £(r p , n). As we will see in the next section, we carry 
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out quantitative analyses not on £(r p , n) directly, but rather two derived quantities: w p (Eq. |j, 
the projection of £(r p , ir) on the r p axis, which does not suffer redshift-space distortions), and 
^vr) (Eq. 0, essentially a cut in £(r p ,ir) at constant r p to measure redshift space distortions). 
We thus compute errors, and the covariance matrix, of the 1-D quantities w p and £(7 t), rather 
than of the 2-D £(r p , n). For each of the samples of Table [I], we create 100 bootstrap realizations, 
and compute w p and £(-7r) for each. Determination of the covariance matrix of errors is then 


straightforward, following [F94a| . [F94a| have discussed the method in detail, showing that it 
gives a good representation of true statistical errors for the correlation function on scales below 
~ 10 /i -1 Mpc, but it tends to overestimate them on larger scales. 

This covariance matrix enables us to fit models for the real space correlation function £(r) 
and the redshift distortions to the observed w p and £(ir) via x 2 (cf., |F94b| ). In practice, the 
effective number of degrees of freedom in the data is smaller than the number of sampled values 
of w p and £(7 t) (i.e., these functions are oversampled), and thus the covariance matrix is singular. 


We therefore follow F94b in using singular value decomposition, which allows the calculation of 
the matrix product in the x 2 function in a robust manner. 


4. REDSHIFT SPACE DISTORTIONS: £(r p , vr) 

4.1. The Observed £(r p , n) 

Figures || @ display the observed £(r p ,7r) both for the complete samples PP-19, PP-19.5, 
PP-20, and for the morphological subsamples E-19.5 and S-19.5. For PP-19, we show the original 
£(r p , 7r), while for all samples, we show a version smoothed by a 3 /i _1 Mpc x 3 h~ l Mpc Gaussian 
to suppress the binning noise and to bring out the global features of the maps. All the statistical 
analyses below are of course carried out on the unfiltered data. 

The contours for PP-19 and PP-19.5 are enormously distorted at very small scales, a 
signature of a high pairwise dispersion <712 (r) at small r, as we shall quantify in § [0. Most 
of this distortion is produced by pairs lying in clusters, in particular in the smallest sample, 
PP-19, which is dominated by half-dozen rich clusters along the Perseus-Pisces chain. The 
small-scale elongations are substantially smaller in the PP-20 sample. This sample has a volume 
of 7.4 x 10 5 /i~ 3 Mpc 3 ; this is quite a bit larger than the volume of PP-19.5, but it does not 
include any more clusters, and therefore the cluster contribution to 1712 ( 1 ) is somewhat diluted. 
This volume is still a factor of 6 smaller than the volume at which Marzke et al. (1995) show - 
using both a COBE-normalized Cold Dark Matter model and a phenomenological model based 
on the observed distribution of Abell cluster velocity dispersions - that ( 712 ( 1 ) stabilizes. For 
samples of volume comparable to those used here, they derive a typical uncertainty on < 712 ( 1 ) of 
~ 180 km s. Although this scatter refers to non-overlapping samples, it does give an explanation 
for the observed difference in a for PP-19 and PP-20. 
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On large scales, the compression of £(r p , 7r) for PP-19 (in particular the £(r p ,7r)= 0.2 and 
£(r p , 7r)= 0 contours) is that expected due to large-scale streaming (cf., |F94b|) . However, the 
Perseus-Pisces supercluster (cf., Figure ^j) lies largely in the plane of the sky, so that the real-space 
distribution of galaxies is intrinsically anisotropic in this sample. Thus it is not a fair sample for 
measuring large-scale streaming. Some part of the observed distortions of the £(r p ,7r) contours 
must be due to infall onto the supercluster itself, as directly observed by Giovanelli et al. (1996), 
who showed that galaxies out to 20 h Mpc from the ridge display infall velocities of the order 
of 1000 km s -1 (see also Willick 1990; Eisenstein et al. 1997). Unfortunately, there is no way to 
disentangle the two effects from £(r p , tt) alone, and we will not discuss the large-scale distortions 
further in this paper. 

Figure |7| shows £(r p , tt) for E-19.5 and S-19.5. The visual difference between the two maps is 
impressive; the ellipticals display a huge small-scale elongation of the contours, while £(r p , t r) for 
the spirals is much more isotropic. This figure demonstrates directly how different the dynamical 
behaviors of the two populations are, and how their real-space correlation functions are mapped 
into redshift space in very different ways. 

We now go on to quantify the real space correlation function, and the redshift distortions, 
from the observed £(r p , it). 


4.2. The Real—Space Correlation Function 


We project £(r p ,7r) onto the r p axis by integrating over the dimension on which the 
redshift-space distortion acts, giving a quantity that is independent of the form and amount of 
the distortion itself, 


w. 


roc roc _ 

<{r P ) = 2/ dy{,{r p ,ir) = 2 dy £ (r 2 + y 2 ) 1/2 
Jo Jo 1 


(5) 


where the second equality follows from the independence of the integral on the redshift-space 
distortions. In the expression on the right-hand side, £ is the real-space correlation function, 
evaluated at r = (r p + y 2 ) 1 ^ 2 . Modelling £(r) as a power law, £(r) = (r/ro) -7 allows us to carry 
out the integral analytically, yielding 


w P {r p ) = r p 


r o 


EdHW) 

m 


( 6 ) 


where F is the Gamma function. We choose ir up , the upper integration limit in Eq. (||), to be large 
enough to give a stable estimate of w p . 


For the PP-19 sample, w p {r p ) is quite insensitive to ir up in the range 20 h 1 Mpc < tt up < 
30 h ^ 1 Mpc for r p < 10 h ^ 1 Mpc. For larger values of r p , w p {r p ) is in fact fairly sensitive to tt up , 







Sample 

Mlim 

dum{h~ l Mpc) 

Morphology 

Ngal 

PP-19 

-19 

79 

All 

1021 

PP-19.5 

-19.5 

100 

All 

852 

PP-19.5-NOP 

-19.5 

100 (No Perseus) 

All 

803 

PP-20 

-20 

126 

All 

577 

S-19 

-19 

79 

Late 

565 

S-19.5 

-19.5 

100 

Late 

481 

S-19.5-NOP 

-19.5 

100 (No Perseus) 

Late 

458 

S-20 

-20 

126 

Late 

333 

E-19.5 

-19.5 

100 

Early 

278 


Table 1: Properties of the volume-limited subsamples analyzed. 


Sample ?’o (h 1 Mpc) 7 


PP-19 

5.95±8;|I 

1 qo+0.04 
1 - yj -0.07 

PP-19.5 

6-95;°i 

1 qq+0.05 
-*--°°-0.07 

PP-19.5-NOP 

« cc+o.34 
d - oo -0.36 

1 qc+ 0.04 
1 -° t, -0.05 

PP-20 

7 05+ 0 ' 47 
' - UO -0.59 

1 79 +O.O 9 

L - ‘ z -0.10 

S-19 

A cc+0.36 
^• oo -0.36 

1 65+ 0 ' 06 
i ’ DO -0.07 

S-19.5 

c cc+0.40 
O-OO_0.45 

1 70 +O.O 7 

L - 1 J -0.08 

S-19.5-NOP 

5 D5 +0 ' 54 
O.UO_o.48 

1 7fi+ 0 - 08 
' D -0.10 

S-20 

5 f)5 +0 ' 61 
°- uo -0.65 

1 qc+0.08 

1,0 —0.09 

E-19.5 

0 qr:+0.75 
5.00-o 76 

2 D5 +0 ' 10 
z - uo -0.08 


Table 2: Best-fit parameters of the real-space correlation function from w p (r p ). 
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but since we are primarily interested in the redshift distortions on small scales, this has no effect 
on our result. 


The observed w p (r p ) and the best-fit power law for the complete samples are shown in 
Figure |], together with likelihood contours on ro and 7 , while the results of the fits are reported 
in Table Note how well the power-law model fits the data|Fj. Error bars are given by the scatter 
over 100 bootstrap realizations and the fit is performed as discussed in § 13.31. There is evidence 


of a growth of the correlation length with increasing sample depth and intrinsic luminosity. This 
is most significant (~ 3<r) between PP-19 and PP-19.5; ro does not grow significantly between 
PP-19.5 and PP-20. This is in qualitative agreement with the results of Iovino et al. (1993) using 
a previous version of this sample, and the results of Loveday et al. (1995) on the APM-Stromlo 
redshift survey, but is in contrast with Hamilton (1988), who found that the luminosity dependence 
of rg was most significant at the highest luminosities. Table |2] also shows a similar trend for the 
spiral-only samples. Thus even within morphological classes, a luminosity dependence of clustering 
does exist (Iovino et al. 1993; cf., their Figure 12). We have checked the sensitivity of these 
results to magnitude errors at the faint end by cutting the PP-19, PP-19.5 and PP-20 samples 
at a corrected magnitude mzw = 15.2, and re-computing £(r p , 7 r), w p (r p ) and the best fit with a 
power-law £(r). This is a fairly conservative selection, reducing each of the three subsamples by 
~ 30% in number (to 882, 740 and 503 galaxies respectively). For these three samples, we obtain 
the following estimates for (ro, 7 ): (5.95^o;go h _1 Mpc, l^l^g®), (6.85^041 h~ l Mpc, 1.90tao<3), 
and (7.45±8;|| hr 1 Mpc, 1.74±g;iJ), respectively. Comparison of these values with those in Table f| 
shows that the results are very robust and that our conclusions are unaffected by any magnitude 
bias affecting the faint end of the Zwicky catalogue. 

Figure ^ shows one of our principal results, the relative clustering strength of early- and 
late-type galaxies, as described by the real-space correlation function. Both the slope and 
correlation length are significantly different in the two samples (Table ||). 

The scale dependence of the relative bias bES of early to late-type galaxies is then simply: 


fcEs(r)= (M v 


(7) 


where 64 is the value at 1 h~ l Mpc and v = ( 7 g — 7 e) 0 ' 5 . We find bEsir) = (2-0 ± 0.4)r ~ 016±0 08 , 
where the error bars have been computed using standard error propagation. At r = 5 /i _1 Mpc, 
we find bES = 1.6 ± 0.4. Hermit et al. (1996) compute a similar relative bias factor from the 
Optical Redshift Survey (Santiago et al. 1995), but in redshift space, finding an average value 
~ 1.5 between 1 and 10/i _ 1 Mpc. Their analysis does not take into account the differences in 
redshift-space distortions between the two classes of galaxies that we have stressed here. Loveday 
et al. (1995) use both the APM catalogue and the sparsely-sampled subsets of galaxies that 


s Note also, however, that w p (r p ) is an integral over £(r), and therefore small deviations from the power-law model 
in the latter function are averaged out in the former. 
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form the APM-Stromlo redshift survey; inverting the angular correlation function w(9) for two 
subsamples limited to bj = 16.57 they find ro = 7.76 ± 0.35 h~ l Mpc, 7 = 1.87 ± 0.07 for early-type 
galaxies and tq = 4.49 ± 0.13 h ^ 1 Mpc, 7 = 1.72 ± 0.05 for late-type galaxies. This is in good 
agreement with our direct estimates from the spatial function. However, the APM-Stromlo data 
are too sparse for being able to compute £(r p , 7 r) directly from the morphological subsamples, so 
that an estimate of the pairwise velocity dispersion 072 (r) cannot be obtained. 

The correlation length we find for all spiral galaxies is significantly larger than that of 
IRAS galaxies in the 1.2 Jy redshift survey, ro(1.2Jy) = 3.76 ^q23 h ~ 1 Mpc, while the slope 
7 ( 1.2 Jy) = 1 . 66^009 is similar. This has the interesting implication that the relative bias of 
spiral galaxies to IRAS -selected galaxies is independent of scale, at least below 10 h~ l Mpc. 
Since IRAS galaxies tend to be of type Sb and later, we have defined a volume-limited sample 
to M = —19, containing 321 galaxies with types between Sb and Irr. In fact, for these we find a 
lower correlation length, ro = 4 . 051 q 75 h -1 Mpc, and a similar logarithmic slope, 7 = I. 55 I 0 I 3 , in 
excellent agreement with IRAS galaxies. 


4.3. The Pairwise Velocity Dispersion 


The quantity £( r p ,7r ) can be expressed as an integral over the product of the real-space 
correlation function £(r), and the distribution function of the line-of-sight components W 3 of 
relative velocities for pairs with separation r, f{w?\r ) ( |F94b ) If y is the component of r along the 
line of sight, then W 3 = Hq{tt — y ) and the integral can be written as (Peebles 1980, F94b) 


/ -(“DO /■ 1 -I -v 

dy{ 1 + f (r 2 + y 2 ) 2 } / [H 0 (tt - y)\r] 

-OO ^ J J 


( 8 ) 


This expression gives a description of the effect of a peculiar velocity field on £(r), but does 
not represent a self-consistent dynamical treatment of the density and velocity fields, which are 
clearly interdependent ( Fisher 1995| ). We do not have any a priori information, therefore, on the 
functional form of the distribution function /. Peebles (1976) first showed that an exponential 
distribution best fits the observed data, a result subsequently confirmed by V-body models (e.g., 
Zurek et al. 1994). With such a choice, Eq. (Isl) becomes 


/ +00 

dy[l + £(r)\ 

-OO 


s/ 2 a 12 (r) 


exp 


-V2Ho 


-K-y 


• 1 V12(r) 
1 H 0 r 


012 (r) 


(9) 


where r 2 = r 2 + y 2 , vi 2 (r) is the mean relative velocity of galaxy pairs with separation r, and 
012 (r) is the pairwise one-dimensional velocity dispersion along the line of sight. 

I I 

F94b show that it is very difficult to model the dependence of vi 2 on the separation r. 
This is made particularly difficult in our case; our sample covers too small a volume to allow a 
determination of the large-scale streaming from the compression of the contours of £(r p ,7r) (cf., 
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0). For this reason, we do not follow [F94b| in a detailed analysis of the mean streaming, and 


instead, we limit ourselves to the simple streaming model introduced by DP83, based on the 
similarity solution of the BBGKY equations: 

F 


VY 2 (r) = -H 0 r 


1 + te)' 


( 10 ) 


We wish to fit Eq. (|J) to the observed £(r p ,7r) to constrain < 712 (r). We are interested in 
particular in < 712 ( 1 ), the pairwise velocity dispersion for scales smaller than 1 h^ 1 Mpc, and thus 
we carry out all fits to the quantity 


£(tt) = f dr p i(r pi Tr) 
Jo 


( 11 ) 


following F94bf J [ In practice, because we have calculated £(r p , ir) in 1 /i _1 Mpc bins, £(tt) is simply 
the value of £(r p ,7r) in the first bin of r p . We assume further that < 712 ( 1 ’) is a weak function of 
separation r ( pP83| ), so that it can be treated as a single free parameter, < 712 ( 1 ). Figure 10 shows 
the results of two-parameter fits of the model of Eqs. (|9|) and ( |l0| ) to £(n) for the E-19.5 and 
S-19.5 subsamples. We use the best-fit values of ro and 7 from Table |] appropriate for each 
subsample; errors and covariances of £(ir) are calculated consistently, as described in § T3. The 
quantity F is very poorly constrained by these data (Figure |l0|): the free-streaming on these small 
scales is quite small. 


We therefore estimate < 712 ( 1 ) for the two cases F = 0 (free expansion of pairs with the Hubble 
flow), and F = 1 (stable clustering). This second case is most probably the one closer to a 
realistic model: Jain (1996) shows that the stable clustering hypothesis (F = 1) should be a good 
approximation at the present epoch for scales of the order of, or smaller than 0.7 h~ l Mpc. We 
thus use the values with F = 1 in our discussion below. 


The value of < 712 ( 1 ) is of order 800 kms -1 for PP-19 and PP-19.5, but drops below 500kms^ 1 
for PP-20, consistent with the more isotropic contours of £(r p ,7r) for this case (Figure |6j). We 
interpret this as due to the smaller effect rich clusters, and in particular the Perseus cluster, have 
on the larger volume of PP-20, as we shall see in the next section. Notice the very significant 
factor of two difference between the < 712 ( 1 ) for early- and late-type galaxies at Mz w = —19.5, a 
dramatic indication of the effect of cluster cores on the determination of < 712 ( 1 ). We now turn to a 
direct demonstration of the sensitivity of < 712 ( 1 ) to the presence of rich clusters in the sample. 


4.4. Stability of <712 for Late—Type Galaxies 

Marzke et al. (1995) have discussed in detail the effect of the contribution of cluster galaxies 
to the small-scale pairwise velocity dispersion. The pairwise velocity dispersion is a pair-weighted 


9 Our definition of £( 7 r) differs from that of F94b| by an unimportant normalization factor - cf., their Eq. (7). 
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statistic, and thus it is heavily weighted in regions of high density, i.e., clusters. Because galaxies 
in clusters have an intrinsically high velocity dispersion, the inclusion or exclusion of clusters 
can have a dramatic effect on Marzke et al. showed that the estimates of fluctuate 

from one sample to the other due to the significant variations in the number of clusters even over 
volumes as large as those of the CfA2 and SSRS2 surveys. Guzzo et al. (1996) showed that < 112 ( 1 ) 
dropped from rs_/ 800 km s 1 to rs_/ 600 km s 1 in PP-19, after removing the Perseus cluster. Thus 
the removal of a single dominant cluster can significantly affect the pairwise velocity dispersion. 


Here we explore further the stability of < 712 ( 1 ), in the case of spiral galaxies. Using spiral-only 
samples, we are in practice filtering out the high-density non-linear regions that would otherwise 
get such high weight in < 712 ( 1 ). The result of excluding the Perseus cluster from the PP-19.5 
sample is visually shown by the changes in £(r p , tt) in the two top panels of Figure |ll|. The 
differences between the contours in the two panels [and the corresponding values of <x 12 (l), 
reported in Table [|, can be compared to those produced by the same operation on the S-19.5 
sample (bottom). While the effect on PP-19.5 is relevant [although less dramatic than it was 
found for PP-19 by Guzzo et al. (1996) for PP-19, due to the larger volume and the consequent 
reduced weight of the Perseus cluster], the two bottom panels of Figure [II] are virtually identical, 
and so are the estimated < 712 ( 1 ) . 


Table [| also gives < 712 ( 1 ) for spiral samples limited to Mz w = —19 and —20. Table ||] indicates 
that < 712 ( 1 ) for spirals lies consistently between 300 and 350 kms -1 . This agrees with the [F94b 
value for IRAS galaxies 317^49 km s 1 , the |Marzke et al. 1995| value for galaxies outside of R > 1 
Abell clusters 295 ± 99kms -1 , and the original determination by DP83| , which undersampled the 
Virgo cluster in CfAl (Somerville, Davis & Primack 1997). 


It is interesting to discuss the similarity of < 712 ( 1 ) for spiral and IRAS galaxies, in the light 
of their different correlation lengths (r 0 — 5.5 h -1 Mpc and ro — 3.8 h~ l Mpc, respectively). IRAS 
galaxies are mostly late-type spirals, and indeed we showed above that if we compute £(r) for this 
subclass, we recover tq ~ 4.0 h~ l Mpc, in agreement with IRAS galaxies (cf., Giovanelli, Haynes 
& Chincarini 1986, and Iovino et al. 1993, who showed that there is a continuity in the clustering 
strength within the spiral class, with Sc’s being less clustered than Sa’s). The similar value of 
< 712 ( 1 ), on the other hand, may simply indicate that the dynamics of IRAS galaxies and all spirals 
are governed by the same fluctuations in the underlying matter density field. Indeed, < 712 ( 1 ) for 
galaxies of type Sb and later is found to be 255^75 5 kms -1 , in statistical agreement with the spiral 
sample as a whole. 


5. SUMMARY AND DISCUSSION 

The main conclusions we have reached in this paper can be summarized as follows. 

• We see very strong small-scale redshift-space distortions in the Perseus-Pisces redshift 
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survey. The distortions are much stronger for early-type galaxies, as one would expect from 
the segregation of morphological types. 


We confirm a mild luminosity dependence of clustering for absolute magnitudes around 
the knee of the luminosity function, Mzw ~ —19.5: the correlation length increases from 
5.95±g-§I h- 1 Mpc for Mz w < —19 to 7.05lg5g/i _1 Mpc for M < —20. These values are 
somewhat higher than the “canonical” value of the correlation length based on the DP83| 
analysis of the CfAl survey, ro = 5.4 h~ x Mpc. This latter value is confirmed by the analysis of 
the ESO Slice Project (ESP) survey (Bartlett et al. 1997), that yields ro = 4.5lg'j ^hr l Mpc, 
while the Las Campanas survey (Lin 1995), gives ro = 5.00 ± 0.14/i _1 Mpc. There are 
probably two reasons for the higher values of ro measured here. First, the “standard” values 
quoted above are estimates of £(r) performed on apparent-magnitude-limited samples. If 
there is, as we have seen, a mild luminosity dependence of clustering, samples that are 
volume-limited at relatively bright absolute magnitudes will systematically measure a higher 
clustering signal. Second, the Perseus-Pisces area is rather rich in clusters of galaxies, so 
that it probably over-emphasizes the cluster contribution to £(r). For comparison, the 
richest cluster in the CfAl volume is the Virgo cluster. The CfA2+SSRS2 sample (Marzke et 
al. 1995) covers part of the PP area and thus has a higher contribution of clusters, resulting 
in ro = 5.97 ± 0.15 /r _1 Mpc. 


• A meaningful comparison of the relative clustering strength of spirals and ellipticals 
can be performed only in real space, i.e. after correcting for the effect of differential 
redshift space distortions. A power-law shape, £(r) = (r/ro) -7 , is a good representation 
of the real-space correlation function between 1 and 10 h~ l Mpc for both ellipticals and 
spirals. Our best-fit estimate of the power-law parameters gives ro = 8.35 !q76 h ^ 1 Mpc, 

7 = 2 . 051 q'q 3 f° r ellipticals, and ro = 5.551q 45 Mpc, 7 = 1.73^QQg for spirals. We model 
the relative bias of the two types of galaxies as a power law with a mild dependence on 
scale, bEs ( r ) = (2-0 ± 0.4) (r/l/i” 1 Mpc) °- 16±0 - 08 i Furthermore, we confirm the continuous 
variation of clustering strength also within the spiral class. For late-type spirals (Sb and 
later), and irregulars, we estimate ro = 4.05io75 Mpc, and 7 = 1.551^3, virtually the 
same correlation function as IRAS galaxies. 

• The quantity 072(1), the measured pairwise velocity dispersion between 0 and l/i _1 Mpc, 
varies considerably between samples of different volumes, going from 855^*5 km s _1 of 
PP-19, to 465l^| kms 1 for PP-20. This variation is consistent with the smallness of the 
volume sampled. 

• The difference in the measured 072(1) between early- and late-type galaxies is remarkable. 
We estimate 0-12(1) = 865^^5 kms -1 for ellipticals and 0-12(1) = 345lg|kms _1 for spirals. 
Contrary to results for the combined sample, the value of 0-12(1) for spiral galaxies alone 
is stable to both changes in the sample volume and the presence of rich clusters. The 
consistency of this value with those measured for non-cluster galaxies (Marzke et al. 1995) 
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and IRAS galaxies ( F94b| ), and the stability among spiral subclasses (for Sb and later types 
we measure 255^5° kms -1 ), suggests that a value <712(1) in the range 300 - 350kms -1 
good estimate of the “temperature” of the galaxy flow outside of virialized structures. 
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Sample 012(1) (F=0) <ti 2 (1) (F=l) 


PP-19 

77r+85 
' ‘ °-65 

855i?i 

PP-19.5 

7Qr+155 

MD -115 

805lJi| 

PP-19.5-NOP 

625l^f 

79^+135 

1 zo -95 

PP-20 

cop:+155 

ozo -115 

4651^1 

S-19 

205 

29511 

S-19.5 

255l| 

345l| 

S-19.5-NOP 

23511 s 

32511 s 

S-20 

415^245 

4851^ 

E-19.5 

815li65 

865l? s ° 


Table 3: Summary of the best estimates of the pairwise velocity dispersion between 0 and 1 h _1 
Mpc, 012(1)) for the two cases F = 0 (free streaming with the Hubble flow), and F = 1 (stable 
clustering). All estimates are in km sec^ 1 . 
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Fig. 1.— The distribution on the sky of all galaxies with Mzwicky < 15.5 after correction for 
extinction. The upper solid line marks the border of the high-extinction region excluded from the 
sample. The large lump of objects near a ~ 3.2 h , <5 ~ 41° is the Perseus cluster. 

Fig. 2.— Galaxy distribution in the volume-limited samples to Mz w = —19.5 and Mz w = —20, 
including all morphological types. 

Fig. 3.— The spatial distribution of early-type and late-type galaxies, volume-limited to Mz w = 
-19.5. 

Fig. 4.— Histograms of the redshift distribution in four representative volume-limited samples. 
Top panels: all galaxies. Bottom panels: early-type and late-type galaxies separately. The dashed 
lines are the distributions expected in the absence of structure. 

Fig. 5.— £(r p , 7r) for PP—19. In this and the following £(r p , n) maps, the heavy contour corresponds 
to £ = 1; for larger values of £, contours are logarithmically spaced, with Alog 10 £ = 0.1; below 
£ = 1, they are linearly spaced with A£ = 0.2 down to £ = 0. The dashed contours represent the 
isotropic correlations expected in the absence of peculiar velocities. The right-hand panel has been 
Gaussian-smoothed with an isotropic filter of width 3 h~ 1 Mpc. 

Fig. 6.— £(r p ,7r) for PP-19.5 and PP-20, Gaussian smoothed. 

Fig. 7.— £(r p ,7r) for the early-type and late-type samples, E-19.5 and S-19.5, respectively, 
Gaussian-smoothed. 

Fig. 8.— The projected correlation function w p (r p ) and results of fits of the power -law model, for 
the three volume-limited samples. The error bars are given by bootstrap resampling. The contours 
give the 68.3%, 95.4% and 99.73% confidence levels on the two parameters taken separately. 

Fig. 9.— As in Figure ||, for early- and late-type galaxies, separately. Note the very significant 
separation in parameter space between the two classes. 

Fig. 10.— Example of full two-parameter fits to £(r p , it) based on the model of Eq. (j|), to show 
how poorly constrained the streaming amplitude F is. Contour levels are as in Figure ||. 

Fig. 11.— The effect on £(r p , n) of removing the largest “Finger of God” from the sample. Right 
panels show the effect of excluding the high-velocity-dispersion Perseus cluster. Top panels are 
for PP-19.5, bottom for S-19.5. Note the change between the two top panels, while for spirals the 
removal of the cluster has very little effect on the measured £(r p ,n). 
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